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ABSTRACT 



C + is a key species in the interstellar medium but its 158 //m fine structure line cannot be observed from ground-based telescopes. 
Current models of fluorine chemistry predict that CF + is the second most important fluorine reservoir, in regions where C + is abundant. 
We detected the / = 1-0 and 7 = 2-1 rotational lines of CF + with high signal-to-noise ratio towards the PDR and dense core 
positions in the Horsehead. Using a rotational diagram analysis, we derive a column density of N(CF + ) = (1.5 - 2.0) x 10 12 cm -2 . 
Because of the simple fluorine chemistry, the CF + column density is proportional to the fluorine abundance. We thus infer the fluorine 
gas-phase abundance to be F/H = (0.6 - 1.5) x 10~ 8 . Photochemical models indicate that CF + is found in the layers where C + is 
abundant. The emission arises in the UV illuminated skin of the nebula, tracing the outermost cloud layers. Indeed, CF + and C + are 
the only species observed to date in the Horsehead with a double peaked line profile caused by kinematics. We therefore propose that 
CF + , which is detectable from the ground, can be used as a proxy of the C + layers. 

Key words. Astrochemistry - ISM: clouds - ISM: molecules - ISM: individual objects: Horsehead nebula - Radio lines: ISM 



1. Introduction 

C + is a key species in the interstellar medium. First, it is an im- 
portant trac er of the neutral ga s where CO has not been able 
to form yet (lLanger et al.ll2010h . Second, it is the dominant gas 
phase reservoir of carbon in the diffuse interstellar medium and 
its fine structure transition CP3/2 - 2 P\/2 157.8 /mi, 1.9 THz) is 
the major cooling mechanism of the diffuse gas. Unfortunately, 
its rest- frame emission cannot be observed from ground-based 
telescopes. With Herschel and SOFIA it is now possible to ob- 
serve this line from space and from the stratosphere respectively, 
but with limited spatial resolution (12" and 15", respectively). It 
is therefore of great interest to find tracers of C + that can be 
observed from the ground at much higher spatial resolution, for 
example with the Atacama Large Millimeter Array (ALMA) or 

the Northern Extended Millimeter Array (NOEM A). 

Th e chemistry of fluorine was studied by iNeufeld et al.l 
(120051) . Although there are still uncertainties in some reaction 
rate coefficients (e.g. CF + photodissociation), the current mod- 
els predict that CF + is present in regions where C + and HF are 
abundant, as it is produced by reactions between these species 
and destroyed mainly by electrons. In these regions CF + is the 
second most important fluorine reservoir, after HF. The ground 
rotational transition of HF, which lies at THz frequencies, was 
first detected with the HIFI/Herschel in the dif fuse interstel- 
lar m edium through absorption meas urements (INeufeld et all 
I2010h and in emission in the Orion Bar (I van der Tak et al ■120121) . 



INeufeld et all (1201 Ol) find a lower limit for the HF abundance of 
6x 10" 9 relative to hydrogen nuclei, providing support to the the- 
oretical prediction that HF is the dominant fluorine reservoir un- 
der a wide variety of interstellar conditions. Unlike HF and C + , 
CF + rotational lines can be observed from the ground. However, 
up to now there has been only o ne detection of CF + towards the 
Orion Bar (INeufeld et al.ll2006l) . 

In this letter we report the detection of CF + at the HCO 
and DCO + peak emission, corresponding respectively to the 
PDR and dense core environmen ts in the Horsehead nebula 
dPetv et al.ll2007l:lGerin et al.ll2009l) . We then infer the gas phase 
fluorine abundance. 

2. Observations and data reduction 

Fig. \T\ displays deep integrations of the J = 1-0 and 7 = 2-1 
low-energy rotational lines of CF + with the IRAM- 30m tele- 
scope centered at the PDR and the dense core, located respec- 
tively at (£RA, SDEC) = (-5", 0") and (20", 22") with respect 
to the projection center, with 49 kHz spectral resolution at both 
frequencies. These observations were obtained as part of the 
Horsehead WHISPER project (Wideband High-resolution Iram- 
30m Surveys at two Positions with Emir Receivers). A presenta- 
tion of the whole survey and the data reduction process will be 
given in Pety et al. 2012, in prep. 

The CF + 7=1-0 and 7 = 2-1 maps displayed in Fig. [2] 
were observed simultaneously during 7 hours of good winter 
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Table 1: Observation parameters for the maps. Their projection center is #2000 = 05M0 m 54.27*, £2000 = -02°28 / 00 // . 



Line 


Frequency 


Instrument 


^eff 


#eff 


Beam 


PA 


Int. Time* 




Noise 


Obs.date 




GHz 








arcsec 




hours 


K(r;) 


K(r mb ) 




HCO lo,i 3/2,2 -Qo,o 1/2,1 


86.670760 


PdBI/C&D 


0.95 


0.78 


6.7 x 4.4 


16 


6.5 


150 


0.09 


2006-2007 


CF + 1 - 


102.587533 


30m/EMIR 


0.94 


0.79 


25.4 





2.5 


88 


0.03^ 


Jan. 2012 


CF + 2 - 1 


205.170520 


30m/EMIR 


0.94 


0.64 


11.4 





3.4 


220 


0.1 8 C 


Jan. 2012 


DCO + 3-2 


216.112582 


30m/HERA 


0.90 


0.52 


11.4 





1.5 


230 


0.10 


Mar. 2006 



a On-source integration time. The noise values are for a spectral resolution of ^0.114 km s 1 and c 0.057 km s l . 



weather (2mm of precipitable water vapor) using the EMIR 
sideband- separation receivers at the IRAM-30m. We used the 
position-switching, on-the-fly observing mode. The off-position 
offsets were (SRA, £Dec) = (-100", 0"), i.e. into the HII region 
ionized by <xOri and free of molecular emission. We observed 
along and perpendicular to the direction of the exciting star in 
zigzags, covering an area of 100" x 100". A description of the 
HCO and DCO + observations and data reductions, which are 
also displayed in Fig. [2l can be found in lGerin et al. and 
Pet v et al] (120071) . Table 1 summarizes the observation parame- 
ters for all these maps. 

The maps were processed with the GILDASQ softwares (iPetyl 

12005b . The IRAM-30m data were first calibrated to the T* scale 
using the chopper- wheel method (iPenzias & Burrusl[l97l . and 
finally converted to main-beam temperatures (T m \>) using the 
forward and main-beam efficiencies (F e ff & 5 e ff) displayed in 
Table 1. The resulting amplitude accuracy is ~ 10%. The re- 
sulting spectra were then baseline-corrected and gridded through 
convolution with a Gaussian to obtain the maps. 

3. Results 

3.1. Line profiles 

Two velocity peaks for the 7 = 1-0 line are clearly seen at both 
positions in Fig.[T] The second velocity peak is marginal for the 
7 = 2-1 line. Table [2]presents the results of dual Gaussian fits. 
The centroid velocity of each peak is significantly shifted be- 
tween the two CF + transitions. We have carefully checked that 
neither peak is a residual line incompletely rejected from the 
image side band. This doubled-peak behavior is unexpected be- 
cause all species without an hyperfine structure detected previ- 
ously at millimeter wavelengths in the Horsehead present a sim- 
ple velocity profile centered close to 10.7 km s" 1 . The only other 
species detected to date with a clear double peak emission profile 
is C + towards the illuminated edge of the cloud (HIFI/Herschel, 
Teyssier et al. 2012 in prep). 

The most obvious explanation would be that the higher ve- 
locity peak corresponds to another line from another species. 
However, there are no other lines in the public line catalogs 
dPickett et al.lll998l: iMiiller et alJl2001b near this frequency be- 
sides CF + . Another possible explanation would be that the two 
peaks correspond to different hyperfine components which are 
caused by the fluorine nuclei. To our knowledge, there are no 
hyperfine structure studies on CF + . However, as the molecule 
is isoelectronic with CO, one can try to rely on 13 CO spec- 
troscopy since the nuclear spin of 13 C is 1/2, as for the fluorine 
nucleus. The magnetic dipolar coupling constant scales approxi- 
mately with the rotational constant and the magnetic moment o f 
the nucleus for l H electronic ground states (iReid & Chull 19741) . 
This allows us to derive a coupling constant of approximately 
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Fig. 1: Deep integrations towards the PDR and core positions. The 
green lines are double gaussian fits. 



110 kHz, and an hyperfine splitting of 165 kHz (110 kHz) for the 
7 = 1-0 (7 = 2-1) transitions, well below the observations. In 
addition, the respective intensities do not follow the theoretical 
predictions. Therefore, this possibility is unlikely. The profile is 
not caused by self- absorption because the CF + opacities are low 
(r ^ 1). We thus attribute the complex line profiles to kinematics 
in the CF + emitting layers. 

3.2. CF + spatial distribution 

Fig.[2]presents the CF + , HCO and DCO + integrated line intensity 
maps. The CF + emission is concentrated towards the edge of the 
Horsehead, delineating the western edge of the DCO + emission. 
A more extended and fainter emission is detected in the CF + 7 = 
1-0 map but not in the 7 = 2-1 map, which has a lower 
signal-to-noise ratio. The intensity peak of the CF + 7 =1-0 line 
coincides with the intensity peak of the HCO emission (shown 
by the green cross ), which traces the far UV illuminated matter 
dGerin et al.ll2009l) . The 7 = 2-1 transition also peaks near the 
HCO emission peak. 



Table 2: Gaussian fit results. 



Line 




Line area 


Velocity 


Width 


T^peak 






Kkms" 1 


km s _1 


km s _1 


K 








PDR 






CF + 7 = 1 


-0 


0.10+0.01 
0.05+0.01 


10.36+0.02 
11.38+0.03 


0.65+0.04 
0.66+0.09 


0.15 
0.07 


CF + 7 = 2 


- 1 


0.25+0.02 
0.04+0.02 


10.62+0.02 
11.39+0.06 


0.67+0.06 
0.38+0.17 


0.36 
0.09 








CORE 






CF + 7 = 1 


-0 


0.03+0.01 
0.01+0.01 


10.47+0.04 
11.32+0.06 


0.60+0.08 
0.45+0.20 


0.05 
0.03 


CF + 7 = 2 


- 1 


< 0.05* 


10.7* 


0.5* 


0.09* 



1 See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the GILD AS softwares. 



* Upper limit for a fixed velocity, linewidth and a r peak = 2cr RMS 
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Fig. 2: Integrated intensity maps of the Horsehead edge. Maps were ro- 
tated by 14° counter-clockwise around the projection center, located at 
(6x, Sy) = (20", 0"), to bring the exciting star direction in the horizontal 
direction and the horizontal zero was set at the PDR edge, delineated by 
the red vertical line. The crosses show the positions of the PDR (green) 
and the dense core (blue). The spatial resolution is plotted in the bottom 
left corner. Values of the contour levels are shown on the color look-up 
table of each image (first contour at 2<x and 2.5cr for CF + 1-0 and 2-1 
respectively). The emission of all lines is integrated between 10.1 and 
11.1 kms" 1 . 



We have checked that beam pick-up contamination from the 
PDR is negligible at the core position (< 7%), even with the 
large beam size (~ 25") of the 30m at 102 GHz. The emission 
then arises in the line of sight towards the core but not necessar- 
ily in the cold gas associated with the core. Indeed, we expect 
the CF + emission to arise in the outer layers of the nebula, delin- 
eating the edge as shown by the maps. This emission is likely to 
arise in the warmer and more diffuse material of the skin layers 
towards the core line of sight. Furthermore, there is a minimum 
of the emission in the CF + 1 - map towards the dense core 
position. This confirms that the CF + emission is associated to 
the diffuse envelope. Molecular emission fro m the lower den- 
sity cl oud surface was alread y mentioned by iGoicoechea et al.l 
(120061) and lGerin et all (120091) to explain the CS and HCO emis- 
sion, respectively. 



^ (4.9 - 6.5) x 10" 10 with respect to H nuclei. Taking the same 
excitation temperature of 10 K, we computed a column density 
towards the core of ^ 4.4 x 10 11 cm" 2 . We consider this as an 
upper limit for the model in Sect. 13. 41 because CF + is found in 
the surface layer which is not taken into account by our unidi- 
mensional photochemical model. 

3.4. CF + chemistry 

CF + is formed through the following chemical path: 

F + H 2 -> HF + H 

HF + C + ^ CF + + H 

and destroyed by d issociative recombination with electrons 
dNeufeld et alJl2006h and by far UV photons: 

CF + +e-^>C + F 

CF + +/*v^>C + +F 

The reactions rates are k x = 7.2 x 10- 9 (r/300)-° 15 cm 3 s _1 
dNeufeld et all 12 005) and k 2 = 5.3 x 10- 8 (r/300)-° 8 cm 3 s" 1 
(iNovotny et al.ll2005l) . The CF + photodissociation rate & p d is not 
known. Nevertheless, assuming a rate of 10" 9 exp(-2.5 Ay) s" 1 , 
we estimate that this contribution is negligible compared to the 
dissociative recombination in low far UV field PDRs like the 
Horsehead. However, it might not be negligible in regions with 
high radiation fields (x - 10 4 - 10 5 ), like the Orion Bar. In the 
following we thus assume that k p a = 0. Because 1) the fluorine 
chemistry is simple, 2) the electron abundance is given by the 
ionization of carbon n(e-) ~ n(C + ) and 3) HF is the major reser- 
voir of fluorine n(HF) ~ n(F), it can be shown that the CF + col- 
umn density is proportional to the fluorine gas phase abundance 
([F]=F/H),/.e. 

N(CF + )- y-\F]nul [cm" 2 ] 

From our CF + observations we find F/H ^ (0.6 - 1.5) x 10~ 8 
in the Horsehead PDR (see Fig. 13 , in good ag reement with 
the solar value (2.6 x 10~ 8 ; lAsplund etaDl2009l) and the on e 
found in the diffuse atomic gas (1.8 x !Q- 8 : ISnow et alJl2007h . 
ISonnentruckeretalJ (1201 Ol) also derived F/H ^ (0.5 - 0.8) x 10" 8 
in diffuse molecular clouds detected in absorption with the 
HIFI/Herschel. 



3.3. Column densities and abundances 

The CF + column density is estimated assuming that the emis- 
sion is optically thin and that the emission fills the 30m beam. 
We infer an excitation temperature of 10 K, based on a rotational 
diagram built with the integrated line intensities of the two de- 
tected transitions. Assuming T ex = 10 K for all rotational levels, 
the beam averaged column density is ^ (1.5 - 2.0) x 10 12 cm -2 
in the PDR. This value is simila r to the column density found in 
the Orion Bar by iNeufeld et al.l (120061) . In the next section, we 
will show t hat the CF + emission arise s at the illuminated edge of 
the nebula. IGoicoechea et al.l (l2009al) found that the [Oi]63 fim 
fine structure line, which also arises at the edge of the nebula, 
was best reproduced with a gas density of nu - 10 4 cm" 3 . Thus, 
assum ing this density and a cloud depth / ~ 0.1 pc (lHabart et al.l 
120051) . the CF + column density translates into an abundance of 
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Fig. 3: Relation between the CF + column density and F/H. 
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In order to understand the CF + abundance profile as a func- 
tion of depth, we used a one-dimensi onal, steady- state pho- 
tochemical model ( Meudon PDR code. iLe Bourlot et all 1201 2l ; 
iLe Petit et all 120061) . The used version of the Meudon PDR 
model includes the Langmuir Hinshelwood and Eley-Rideal 
mechanisms to describe the formation of H2 on grain surfaces. 
The physical conditions in the Horsehead have already been con- 
strained by our previous observational studies and we keep the 
same assumptions for the steep density gradient (displayed in 
the upper pan el of Fig. [4]), radiation field (x = 60 times the 
iDraind (Il978l) mean interstellar r adiation field), elemental gas- 
phase abundances (see Table 6 in iGoicoechea et al1l200 9b) and 
cosmic ray primary ionization rate (£ = 5 x 10" 17 s" 1 per H2 
molecule). We used the Ohio State University (osu) pure gas- 
phase chemical network, and included fluorine adsorption and 
desorption on grains. 

The predicted CF + , HF and C + abundance profiles are shown 
in Fig. |4] (b). HF and CF + abundances decreases rapidly for 
Ay > 1. The model is in good agreement with the observed CF + 
abundance in the PDR, shown by the horizontal bars. The model 
predicts that there is a significant overlap between CF + and C + . 
Moreover, the abundance ratio between these two species re- 
mains quite constant along the illuminated side of the cloud, i.e. 
A v < 4, as shown in Fig.|4](c). The CF + emission arises in the 
outermost layers of the cloud (Ay ~ 0.5), which are directly ex- 
posed to the far UV field and where the gas is less dense. The 
predicted spatial distribution of the CF + emission is shown in 
Fig. |4] (d). We expect a narrow filament (~ 5") shifted in the 
illuminated part of the PDR with respect to the HCO emission, 
which h as already shown t o trace the far UV illuminated molec- 
ular gas (iGerin et al .1120091) . 

4. Conclusions 

We have detected the / = 1-0 and J = 2 - I rotational lines 
of CF + with high signal-to-noise ratio towards the PDR and core 
positions in the Horsehead. We have also mapped the region, and 
find that the emission arises mostly at the illuminated edge of 
the nebula (PDR), but it is also detected towards the dense core 
arising from its lower density skin. CF + is unique as its column 
density is proportional to the elemental abundance of fluorine. 
In the Horsehead PDR we find N - (1.5 - 2.0) x 10 12 cm" 2 
and infer F/H (0.6 - 1.5) x 10" 8 . Our model of the fluorine 
chemistry predicts that CF + accounts for 4-8% of all fluorine. 
CF + is found in the layers where C + is abundant as it is formed 
by reactions between HF and C + . In these reg ions the ionization 
fraction is high (see IGoicoechea et al.ll2009bl) and CF + destruc- 
tion is dominated by dissociative recombination with electrons. 
The CF + emission has two velocity components. The possibil- 
ity that we are resolving the hyperfine structure is unlikely but 
corresponding theoretical or experimental study would allow to 
derive the velocity structure unambiguously. Although the CF + 
line profile is not exactly the same as the C + line profile, they are 
the only species in the Horsehead with a double-peaked profile 
of kinematic origin measured to date. The complex line profile 
of both CF + and C + therefore confirms that they trace the gas 
directly exposed to the far UV radiation, which shows a com- 
pletely different kinematics than the following layers traced by 
other species, like HCO. We therefore propose that CF + can be 
used as a proxy of the C + layers which can be observed from the 
ground. We will check this by comparing with a HIFI/Herschel 
map of the C + emission in the Horsehead. 
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